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ABSTRACT 

Panchromatic spectral energy distribution (SED) from the ultraviolet (UV), 
optical to infrared (IR) photometry of NGC 628, combined with the evolutionary 
stellar population synthesis, is used to derive the spatially resolved age, metal- 
licity and reddening maps. These parameter distributions show that the bulge 
of this galaxy is a disk-like pseudobulge, which has the Sersic index close to the 
exponential law, rich gas, and a young circumnuclear ring structure. We also dis- 
cover the disk has two distinct regions with different radial age and metallicity 
gradients. The inner region is older and has a much steeper age gradient than 
the outer region of the disk. Both these two regions and the central young struc- 
ture can be seen in the radial profile of the optical color. Based on the age and 
reddening distributions, we consider that the pseudobulge and disk are likely to 
have grown via the secular evolution, which is the redistribution of mass and en- 
ergy through the angular momentum transport caused by the non-axisymmetric 
potential of the spirals. However, possible gas accretion events could affect the 
outer region of the disk, due to abundant Hi gas accumulating in the outer disk. 

Subject headings: galaxies: evolution — galaxies: individual (NGC 628) — galax- 
ies: photometry — galaxies: stellar content 



1. INTRODUCTION 

Galactic evolution at the early universe was dominated by the hierarchical clustering 
and merging that are violent and rapid, secular evolution will be most important in the 



far future, and now both processes might be equally fundamental ( iKormendy fc Kennicutt 
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20041 ) ■ Recently, more and more observations have revealed that the bulge with the Sersic 
inde x approximating the exponeiitial law might be grown via the secular evolution of the 



disk (iKormendv 



Ganda et al. 



200 



1993 

T 



Boker et al 



20021 : iKormendv fc Kennicuttll2004j : iFisher fc Drorvll2008 



Fisher et al.l 120091). This kind of b ulge is called pseudobulge, which is 



very common in late-ty pe galaxies (iBoker et al.l |2002| ) . Pseudobulges have smaller Sersic 
indices than that of the de VaucouleursI ( 1948 ) law in classical bu lges, central gas concen 



t rat io n, active star formation, nuclear rings, bars and/or spirals (IKormendy fc Kennicutt 



20041 ). Non-axisymmetric gravitational potential by the bar, oval, and/or spiral structures 
rearranges the angular momentum and mass, makes gas infall, and triggers the central sta r 
formation to create a pseudobulge (IRegan &: TeubenI l2003l : iKormendy fc KennicuttI 120041 ). 
The gas infall from the disk and the formation of the bulge are secular. In addition to the 
secular evolution, distant encounters, gravitational interactions in clusters, and gas accre- 
tions from the neighboring galaxies can affect the formation and evolution of the bulge an d 
disk dToomre fc Toomrelll972l : iNoguchilEgSsI : ICortes et allbooei iBournaud fc Combeslliooi ) . 
Thus, isolated nearby galaxies, which have large apparent sizes and are far from the tidal 
interactions of other galaxies, are excellent objects to investigate the secular evolution of the 
disk and formation of the bulge. 

Spatially resolved age, abundance and reddening distributions can help us to understand 
the galactic evolution. They provide some basic information about the star format i on his tory. 



chemical comp ositions and evolution, and interstellar medium. Both iKong et al.l ( l2000l ) and 



Li et al.l ( 120041 ) adopted different evolutionary population synthesis (EPS) models to analyze 
the structure and evolution of a nearby Sab spiral galaxy (M81) using the photometric data of 
multiple intermediate-band filters from 3000 A and 10000 A in the Beijing- Arizona- Taiwan- 



to study the properties of stellar populations in galaxies ( 


Fioc & Rocca-Volmeranee 


1997; 


Leitherer et al. 


1999; 


Bruzual & Chariot 


2003; 


Kotulla et al. 


2009) and it has been applied 


in many fields 


such as star clusters, galaxies, and galactic clusters ( 


Abraham et al. 


1996; 



Kong et al.ll2000l : Ma et al.ll2009af ). 



In this paper, we will use the EPS model to derive the population properties (age, 
metallicity and intrinsic reddening) and gain the clues for the formation and evolution of a 
nearby face-on SA(s)c galaxy, NGC 628. This galaxy, al so named M74, has an apparent size 
of about 10' and it is isolated (IKamphuis fc Briggslll992l ). Light decompos ition and color dis- 



tributions reveal that this galaxy has a disk-like bulge (i.e. , pseudobulge) ( iNatali et al.lll992 



Cornett et al. 


1994; 


Sanchez et al. 


2011; 



band observations taken with the 60 /90 cm Schmidt te lesco pe of National Astronomical Ob 



servatories of China as also used by iKong et al.l (120001 ) and 



Li et al.l (120041 ). there are plenty 



of on-line archival data for this galaxy from various surveys and telescopes, such as the 
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Galaxy Evolution Explorer (GALEX), Two Micron All Sky Survey (2MASS), and Spitzer 
space telescope. These panchromatic photometric data-set from UV to IR can considerably 
aid us to obtain more accurate estimations of the properties (age, metallicity and intrinsic 
reddening) . 

Previous studies on NGC 628 present various properties in respect of colors, stellar 
populations, abundance, and reddening, etc., which can be compared with our results. The 
color of the whole gala xy is very blue, rang in g from the blue s t of M 33's colors to the bluest 
of M81's in UV - 1/ JCornett et al.lll994h . ICornett et aP (|l994h reported that the disk 
might have undergone significant star formations over the past 500 Myr by co mparing the 



UV /o ptical colors of stellar population synthesis models with those of NGC 628. iNatali et al. 



(I1992I ) highlig hted that two differe nt stellar populations dominate the inner and outer disk 
of the galaxy. ISanchez et al.l (1201 ll ) presented the age and abundance profiles calculated by 
the integral field spectroscopy in an unprecedented high spatial resolution. They found a 
young circumnuclear star formation region near the galactic center. Distinct comp onents 
with d iffere nt ages and metallic i ties c an be easily discriminated in those profiles. [Hodge 
(119761 ) and iKennicutt fc Hodgd (Il980l ) both used Ha images to obtain several hundreds 
of Hll regions. Numerous Hll regions, tending to lie along the spiral arms, indica t e that 
a n umber of stars have been forming recently in this nearby galaxy. Both iTalentI (119831 ) 
and iBelley &: RoyI (Il992[ ) discovered a notable radial metallicity gradient ac cording to the 
spectr ophotometry and narrow-band imaging of Hll regions. Reddening in iBelley &: Roy 
(119921 ) does not disp l ay an y trend as a function of the radial distance from the galactic 



center. iRegan et al.l ( 120061 ) investigated the 8 /im emission from the polycyclic aromatic 
hydrocarbons (PAHs) for some disk galaxies (including NGC 628) and declared the 8 yum 
PAH surface brightness could be used as an approximate tracer of the interstellar medium. 

The paper is organized as follows. In Section m the observations from the BATC survey 
and the archival data are described. Detailed data reduction is given in Section El The 
introduction of the EPS model and fitting method is provided in Section |H Decomposition 
of the surface brightness profile and the distributions of age, metallicity and reddening as 
well as relevant comparisons are presented in Section [51 Discussions and conclusions are 
made in Section [HI and [3 respectively. 



2. OBSERVATIONS 



NGC 628, as one of the nearby galaxies in the BATC sky survey ( iBurstein et al.[[l994[ ). 
was observed by the 60/90 cm Schmidt telescope, which is deployed at the Xinglong Station 
belonging to National Astronomical Observatories of China. A 2048x2048 Ford Aerospace 
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CCD with the pixel scale of 1.7" (15 fim pixel"^) is mounted at the focal plane of the 
telescope (the focal ratio is f/3). The size of field of view (FOV) is about 58' x 58'. For the 
photometric system, 15 intermediate-band filters, covering the wavelength range of 3300 - 
10000 A with bandwidths of about 200 - 300 A, are well designed to avoid strong emission 
lines of the sky light (IFan et al.l 119961 ) . The observation mode is divided into two parts: 
deep exposures for observations of the objects; short expo sures at photome tric nights for 
both the objects and standard stars to do flux calibrations (jZhou et al.l 120011 ) . The normal 
observations for all filters started in 1995 November and ended in 2004 February. Because of 
the bad image qualities of a, b, and c filters, we reobserved NGC 628 in 2009 winter for these 
three bands. Table [1] summarizes the basic parameters of filters, the total exposure time for 
each band, and some information of calibrations. 

Figure [T] displays the combined images of NGC 628 in three BATC bands with wave- 
lengths centered on 3890, 6660, and 9190 A. The galaxy is a grand-designed spiral galaxy 
with two symmetric inner arm s and multiple long and c ontinu ous outer arms (classified as 
Type 9 in the arm classes of lElmegreen &: Elmegreeru (119871 )). We can see in those fig- 
ures that two clear spiral arms extend from the center to the disk. Some Hll regions can 
be clearly recognized near the spiral arms in the 6660 A band (this band covers the Ha 
emission line) and in the 3890 A (near the ultraviolet regions where massive young stellar 
populations generate strong emissions), while the stellar mass distribution of the galaxy is 
visi bly reflected in the near-infrared wavelength of 9190 A. The disk inclination is about 
6° (IShostak fc van der Kruitlll984j : iKamphuis fc Briggslll9921) and the po sition angle of the 
major axis is about 25°. The redshift is about 0.0022 (iHuchra et al.lll999l ) and the redshift- 
independent distance is about 8.6 Mpc ( Herrmann et al.l 20081 ). This distance together with 
the disc inclination and position angle is adopted throughout our paper. 

In addition to the optical data of the BATC survey, the archival data of the ultraviolet 
and infrared observations are collected to obtained the panchromatic spectral energy dis- 
tribution of the galaxy from the on-line database. In the Nearby Galaxy Survey (NGS) of 
GALEX, NGC 628 was imaged in the far-UV (FUV) and near-UV (NUV) bands, providing 
substantially accurate ultraviolet observations wit h photometric prec i sion up to 0.05 and 0.0 3 
mag in the FUV and NUV bands, respectively ( iMartin et al.ll2005l : iMorrissey et al.l 120071 ). 
Due to its bad fitting to the model of stellar population synthesis, whose spectra l libraries 



l ack ex treme horizontal-branch stars in the old stellar populations as mentioned in iMa et al. 



( l2009bl ). the GALEX FUV band is excluded in the subsequent analysis. XMM Optical Mon- 
itor (XMM Q M) also observed t his galaxy in three other ultraviolet bands (UVWl, UVM2, 
and UVW2) jKuntz et al.ll2008h . The signal to noise ratio (SNR) of the UVW2 band is 
extremely low, so this band is not used in our study. The UVM2 band together with UVWl 
reinforces the observation data in the near-UV region. Near-infrared observations of NGC 
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628 in 2MASS J, and Ks bands (ISkrutskie et al.ll2006l ) are gained. Since only a small 
core in the band (radius less than half an arcmin) has enough S NR and the public data 
of this band from the 2.3 m Bok telescope (iMcCarthy et al.ll200ll ) is available, we replace 
the 2MASS Ks observation by th e image of the B ok telescope for a detailed descrip- 
tion of the data reduction refer to iKnapen et al.l (12004r)'). Infrared o bservations in 3.6 and 
4.5 /um of the Infrared Array Camera (IRAC, see iFazio et al.l (120041 )) of the Spitzer space 
telescope (jWerner et al.ll2004j ) expand the infrared data of NGC 628. Other bands in this 
telescope, such as 5.8, 8.0, 24, 7 0, and 160 /im, are located in the radiation regions of gas 
and dust (IKennicutt et al.ll2003l ). Our study puts emphasis on the star light emitted from 
the stellar content, so these bands are not used in the analysis of the stellar populations. 
The ultraviolet, optical and infrared information of those telescopes and surveys including 
the instruments, images, and corresponding references is listed in Table [2J 



3. Data Reduction 
3.1. Processing Flow of BATC Images 

Subtraction of the bias and dark current and the f l at-fiel d correction were performed 



by the PIPLINE 1 procedure as described in iFan et al.l (119961 ) once the CCD frames were 



observed. Astrometric information was appended to the image header during this step. Deep 
frames of each band were combined after shift and rotation to the same image center and ori- 
entation according to the star positions. Cosmic ray hits and bad pixel effects were fixed up 
during the image combination. Short exposures for both the object and the standard stars 
were taken at photometric nights. We use these images to perform the flux c alibration for 



the a bove deep combined image as follows: 1) Four Oke-Gunn standard stars (jOke fc Gunn 



19831 ) are extracted in short exposure images of the standard stars and t heir instrument al 



magnitudes are measured by a stellar photometric procedure (DAOPHOT; IStetsoru (119871 )): 
2) Comparing the instrumental magnitudes and calibrated AB magnitudes of these standard 
stars, we get atmospheric extinction coefficient and instrument zeropoint; 3) Tens or hun- 
dreds of unsaturated bright stars are selected as secondary standard stars in a short exposure 
image of the object, which is observed under the same photometric condition as the standard 
stars; 4) according to the extinction coefficient and zeropoint, the instrumental magnitudes 
of these stars are converted to fiux-calibrated AB magnitudes. 5) The average difference 
between the instrumental magnitudes of the same stars in the deep combined image and 
their AB magnitudes are calculated; 6) with this difference, an analog-to-digital unit (ADU) 
can be converted calibrated flux density unit of ergs s~^cm~^Hz~^. Calibration error of each 
short exposure for the object contains two items: one is from the magnitude estimations of 
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the standard stars and the other is from the differences of those secondary standard stars 
between instrumental magnitudes and AB ones. The total calibration error for each band as 
shown in the last column of Tabled] is the root mean square (RMS) of the calibration errors 
in all short exposures, the number of which is presented in the seventh column of Table [T] 



3.2. Background Subtraction and Image Scaling 

Sky backgrounds should be removed from the mosaics of all bands in order to acquire the 
intrinsic flux of the galaxy. For the GALEX mosaics, background images are provided, but 
the background near the center is overestimated if we check their sky background intensity 
images. The mosaics from the 2MASS, Bok Ks band, and Spitzer were constructed after 
removal of the sky background by their own techniques (see references in Table [2]). We obtain 
the smoothing sky background maps for all other mosaics by the polynomial fitting method, 
based on th e remaining background pixels after masking the source signals extracted by 



SExtractor (IBertin fc Arnoutslll996l ) and the galaxy by circling it with temperate apertures 



for different bands. These fitted sky backgrounds are then removed from the original mosaics. 

Due to the diverse pixel scales ranging from 0.75" to 1.7 " and different directions of the 
mosaics, we adjust them to the same scale and direction as the BATC images (i.e., 1.7" per 
pixel, top of north and left of east). Take the Spitzer mosaic for example, the pixel size is 
0.75" but that of the BATC mosaic is 1.7". We first calculate the equatorial coordinate for 
each pixel of the BATC image, and then find the corresponding pixel position in the Spitzer 
mosaic according to its astrometry. After counting how many Spitzer pixels (can be decimals) 
are covered by the area of 1.7x1.7 arcsec^ centered at this pixel position, we sum all the 
fluxes of those pixels and subsequently the adjustments of both the pixel scale and direction 
are completed simultaneously. During the adjustments, we transform the pixel values of all 
bands to AB top-of-atmosphere fluxes in unit of 10~^°ergs s~^cm" "^Hz-^ For the Bok 2.2 
//m image from the NASA/IPAC Extragalactic Database (NED), since its astrometry is not 
the same as ours and the magnitude zeropoint is uncertain, we do something special: several 
stars with astrometric positions in the BATC image are used to update the astrometry of 
this image; mosaic fluxes are recalibrated to the AB magnitude system by some selected 
bright star in the same fleld from the 2MASS catalogues. 

Figure H] displays the true color image of NGC 628, combined by three scaled mosaics 
of the NUV, narrow Ha (from NED), and 8 //m bands, which are chosen to present the 
distributions of young stellar population and dust. From this image, we can flnd that the 
whole disk is very blue, implying much younger populations in its disk than in the bulge. 
Numerous Hll regions, identifled as luminous spots, mostly lie along the spiral arms. Light in 
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Hll regions comes from both the old galactic background stellar population and young stars, 
such as O and B type, which have been born recently as a simple stellar population. These 
regions are more complex than the galaxy itself and they will be specially treated in our 
series papers. So pi xels belonging to 3 76 Hll regions in all mosaics are obliterated according 
to the catalogue of iFathi et al.l (120071 ) . We also find numerous filaments and substructures 
generated by the PAH emission in red color of Figure El indicating the gas/dust distribution 
across the disk in addition to the strong PAH radiation along the spiral arms. 



3.3. Image Smoothing and SED Extraction 

Mosaics are smoothed for reasons that: 1) different astrometric precisions for telescopes 
lead a single pixel to represent adjacent slightly different parts of the galaxy; 2) transforma- 
tions and rotations induce position errors; 3) different observational conditions cause diverse 
image qualities (e.g., CCD noises and seeing); 4) SNR varies from the center to the edge of 
the galaxy so that the photometric accuracy decreases in low SNR areas. Smoothing can 
improve SNR but in return reduce the spatial resolution of the galaxy. We use the boxcar 
averaging method to smooth all the mosaics as expressed in the following formula: 

^ /+f J+f 

K=I-^ L=J-f 

where / and J specify the pixel position in the image, M/ j is the smoothed flux of this pixel, 
Fk^l is the original flux of any given pixel {K,L), and w is the smoothing window width 
(so called boxcar width), depending on the SNR of the specified pixel in the BATC j band 
image. The SNR {R) of a pixel is defined as 

p^_S _ S 

where S is the pixel flux, N is the noise, and 6 is the standard deviation of the global sky 
background. Here, SNR is similar to the nominal definition in a CCD image, and we only 
use it to calculate the boxcar width and estimate relative photometric errors for different 
bands. The window size for each pixel is determined hy w = 2mm{[^], 5} + 1 where [^] 
is the minimum integer larger than ^ and Rm is set to 10. The formula promises w to be 
an odd number to obviate the smoothing bias due to asymmetry. It means that if SNR of 
the pixel is equal or larger than 10, the minimum smooth width is 3 pixels and if SNR is 
equal or less than 1, the maximum width is 11. 

After smoothing the mosaics, we extract the observed SED for each pixel. The observed 
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SED includes AB magnitudes of 23 bands from UV to IR, which are calculated as 

= -2.51ogio^ - 48.6, (1) 
as long as the flux signal S is larger than 36. The magnitude error (a) is determined by 

a = 2.51ogio(^^^) = 2.51ogio(l + j). (2) 

At least 15 out of 23 bands with valid magnitudes are required for each SED to fit with 
model ones as explained in the next section. As a result, we get 33242 SEDs and their errors 
for the whole galaxy (see Table [3] about the SEDs of 6 randomly sampled pixels). 



4. STELLAR POPULATION SYNTHESIS AND FITTING METHOD 



Evolutionary population synthesis has become a popular technique in studying dif- 
ferent properties of stellar populations and the evolution histories o f galaxies, since there 



were some significant progresses on the theory of stel l ar evolution ( Renzini &: VoU 



Maeder fc Mevnetl Il989l : Ixippenhahn fc WeigertI llQQol: ISchaller et al.l llOoS 



especially, the diverse phases of evolution (llben Sz Renzinil 



Vassiliadis Sz Wood 1993; Groenewegen fc de JonJl993h and the completeness of stellar ob- 



serva tional and theoretical spectral libraries ( iLe Borgne et al. 



198^ 



1981 



Bertelli et al. 



Schoenberne 



1983 



19971 ). Most f ashionable synthesis m odels, like PE GASE fjFioc fc Rocca-Volmerangel 119971 ). 



StarBurst99 flLeitherer et all Il999l ). GALEXEV flBruzual fc Charlo 



2003 



Pickles 



1998 



Leieune et al. 



jKotuUa et a l.l l2009l). are used to study globular clusters (iMa et al.l 



(Abraham et al. 



20001). galaxies ([Kong et al.ll2000l : IPapovich et al.ll200ll : iLi et al.l 120041 ) . and galactic clusters 



200 



i 



and GALEV 



2009al : iBarmbv et al. 



19961 ). etc. 



In the series papers on the study of nearby gal axies based o n the photometric obser- 
vations of the BATC multicolor sky survey program, iKong et al.l (|2000l ) first quantified the 
chemical a bundance, age and red d ening di stributions of M 81 using the synthesis models of 



SSP (from iBruzual A. &: Chariot! (119931 )). iLi et al.l ( 120041 ) then used those spectral energy 



distributions (13 intermediate band filters) to gain the age and metallicity distributions of 
M81 with the help of PEGASE, which contains exponentially decreasing star formation rates 
(SFRs). Comparison of the above two models are described detailedly in their papers. We 
intend to apply the PEGASE model to investigate the distributions of the stellar population 
properties in NGC 628. 
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4.1. PEGASE Model 



PEGASE is a spectrophotometric evolution model for starburst and evolved galaxies of 
the Hubble sequence attributed to the extension to the NIR atlas of synthetic spectra and 
revised stellar libraries including cold stellar parameters, stellar tracks of asymptotic giant 
branch (AGB), and post-AGB phase (iFioc fc Rocca-Volmerangd 119971 ) . The stellar tracks 
they adopt are mainly from the "Padova" group, which contain a wide range of chemical 
abundances, Z = 0.0001, 0.0004, 0.004, 0.008, 0.02, 0.05, and 0.1 with Y = 2.5Z + 0.23 
{Zq = 0.02, where Z is the metallicity abundance and Y is the helium abundance). The 
initial stellar mass in the tracks ranges from 0.6Mq to I2OM0, where in the track of Z = 
0.1 pseudo-tracks with mass larger than 9Mq are generated from the corresponding masses 
in the Z = 0.02 and Z = 0.05 tracks. 

The stellar spectral libr aries for optical wavelen gth (3130-10800A) come from the stellar 
spectrophotometric atlas of iGunn fc Strykerl (119831 ) which includes 175 stars with complete 
ranges of spectral type and luminosity class. For the far-UV (1230-3200A), stellar spec- 
tra are extracted from the International Ultraviolet Explorer (lUE) ESA/NASA libraries 
( Heck et al.l [l98J). In the extreme- UV (220-1230 A), s pectra f o r the effective temper- 
a te Teff < 50000 K are ca lculated using the models of iKuruczl (119921 ). The models of 



Clegg fc Middlemasd (119871 ) are used at all wavelength for hotter stars (Teg > 50000 K). 
For cold stars dominating the NIR (1-5 fim) of the spectra, observational libraries when 
possible and synthetic spectra otherwise are adopted. Blackbody radiations supplement the 
N IR wavel e ngth range of hotter stars. In the mid- infrared (> 5yum) the analytic extension 
of lEngelkd (119921 ) is used for stars colder than 6000 K. 



Simple stellar population (SSP) is regarded as a simple system whose stars are born at 
the same time and with the same initial chemical compositions (typical example: globular 
clusters). SSPs can be calculated in PEGASE with initial mass functions (IMFs), stellar 
libraries , and stel l ar evo lutionary tracks. PEGASE provides several commonly applied IMFs, 
and the ISalpeten (119551 ) law with stellar mass larger than 0.1 Mq and less than 120 Mq is 
used to generate SSP models in this paper. 

Galaxies are much more complex systems than globular clusters. Evolutionary com- 
posite stellar population (CSP) can more exactly describe the evolution history of galaxies. 
CSP is considered as a superimposition of SSPs of different ages. The integrated spectrum 
of CSP is synthesized with an IMF and a SFR for every metallicity. We assume each pixel 
of the mosaics to be such a CSP system that is composed of stars formed in different peri- 
ods. There are three kinds of SFRs in PEGASE and the exponentially decr easing SFR with 
_15Gyr is used in our paper. Here, r = 15Gyr is typical for Sc galaxies (IBolzonella et al. 



r 



2OOOI). 
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Another feature of PEGASE is that the model takes into account the nebular emission 
(continuum and lines) generated by the ionized gas in the star-forming regions. Nevertheless, 
the regimes used to synthesize the spectral energy distribution for Hll regions are very 
outdated. We do not consider the nebular emission in the model spectra. 



4.2. Fitting Method 



For model spectra, both age from to 20 Gyr and metallicity from 0.0001 t o 0.1 are 
i nterp olated by the PEGASE code in the logarithmic space with steps of 0.01. The ISalpeter 
( 119551 ) IMF and exponentially decreasing SFR are used to convolve t hese spectra to create 
CSPs of the specified metallicity. We use the dust extinction model of ICardelli et al.l (119891 ) 
to redden the CSP spectra with a reddening step of 0.01 ranging from to 1.0. By convolving 
the reddened CSP spectra with the filter transmission curves of all bands, we can obtain 
the model SEDs set in the ranges of different ages, abundances, and reddenings. The model 
SED contains convolved AB magnitudes of 23 bands computed as 



m 



csp 



-2.51og 



10" 



J^F^{t,Z, E)T,{X)dX 



-48.6, 



(3) 



where i is the specified band index, Tj is the corresponding filter transmission curve, F{t, Z, E) 
is the CSP spectrum of the specified age t, metallicity Z and reddening value E in E{B — V), 
and m^^^ is the resultant synthetical AB magnitude. By comparing the observed SED (see 
Section 13. 3p with the synthetical model SEDs, we can simultaneously derive those three pa- 



Kong et al. 


(2000 


) and 


Ma et al. 


(2009b) 



23 



x\t,Z,E) = Y^ 



\m. 



obs 



— m 



csp 



{t,Z,E)f 



where m°^^ is the observed magnitude of the specified filter in the observed SED as shown 
in Equation m (Tj is the observed magnitude errors as given in Equation |51 and rn^^{t, Z, E) 
is the corresponding synthetical magnitude in the model SED as shown Equation [31 

The degenerate effect of age, metallicity, and dust reddening exists when we match 
the observed SEDs with the model ones predicted by EPS, since the ways of these three 
parameters to affecting spectra are similar. However, dust and gas absorb the UV radiation 
of star light and reradiate out in the infrared region. The multi-band photometric data from 
ultraviolet to infrared may aid in degrading this degeneracy. Before being fitted by the above 
method, the observed SEDs are corrected with the foreground reddening of the Galaxy by 
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using the same reddening law of ICardelli et al.l ( 119891 ) . NGC 628 is located far away from 
the Galactic disk as the Galactic latitude of its center is b = -45.7° (R.A.: 24.17° and 
Dec: 15.78°). So th e foreground Galactic reddening is small, about 0.07 mag in E{B — V) 
JSchlegel et al.lll998h . 



4.3. Fitting Results of Several Selected Pixels 

We randomly sample 6 pixels from different parts of NGC 628 to illustrate the SEDs, 
fitted parameters, and the fitting goodness: No. 1 near the galactic core. No. 2 from the 
inner region of the disk. No. 3 from the outer region of the disk. No. 4 from the inner spiral 
arm. No. 5 from the outer spiral arm, and No. 6 close to an Hll region (see Table [3]). The 
positions of these six pixels are marked in the true color map of Figure [2] and the observed 
SEDs, matched model SEDs, and the corresponding model spectra are plotted in Figure 
lU Excellent match to the model spectrum can be seen in this figure. We might coarsely 
conclude that: the core is oldest; age becomes younger distant from the center; spiral arms 
are younger than the rest components except for the places near the Hll regions with much 
younger ages and richer abundances. 



5. Structure and Stellar Population Analysis of NGC 628 

Based on the spectral synthesis model and fitting method as described in the previous 
section, we obtain spatially resolved age, metallicity and intrinsic reddening distributions of 
NGC 628 after fitting the observed SED of each pixel with the models. In order to analyze 
the properties of different components, we first decompose the galaxy into two components in 
both the optical and near-infrared surface photometry. Then the stellar population analysis 
of the whole galaxy and its components is presented in the following sections. 



5.1. Bulge-disk Decomposition of the Surface Brightness Profiles 



As generally known in spiral galaxies, the bulge and disk exhibit many physical and 
dynamical differences. The bulge is usually brighter, older, and dynamically hotter than 
the disk. In order to analyze different components of NGC 628 with respect to the two- 
dimensional features, we disassemble the galaxy into two main components. These two 
components are obtaine d by fitting the azimuth ally averaged surface brightness with the 
exponential disk profile (IFreemanl Il970l ) and the ISersid (119681 ) law following the iterative 
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decomposition procedure of iKormendyl (Il977l ). In this decomposition method, two separate 
regions of the disk and bulge are chosen in the beginning. Then, a least square fitting of 
an exponential law to the radial surface brightness in the disk range is performed. The 
calculated exponential disk contribution is extrapolated to the bulge range and subtracted 
from t he observed p rofile to get a first estimate of the bulge c ompo nent, which has been fitted 
by the ISersid ( 1l968l ) law. In the same way, this fitted ISersid (119681 ) law is extrapolated to the 
disk range and subtracted from the observed profile to get an estimate of the underlying disk, 
which be fitted by the exponential law again. The procedure repeats until all the parameters 
of two laws converge within a specified accuracy (10"'^). 

The radial profile is the corrected azimuthally-averaged values, given the disk inclination 
of 6°, the position angle of the major axis of 25°, and the distance of 8.6 Mpc as previously 
mentioned. The value at each radius is defined as the average of the pixel values (here, 
surface brightness) within a suitably specified annulus of that radius centered on the galaxy 
nucleus. The profile error is calculated as the standard deviation of those pixel values. All 
radial profiles of different parameters (age, color, metallicity, and reddening) in the rest of the 
paper are computed in the same way with those adopted parameters (inclination, position 
angle and distance). 

Left panel of Figure [3] shows the radial surface brightness profile of the BATC d band 
(filled squares with error bars), whose effective wavelength is close to the broad B band, 
and the decomposed bulge and disk components as shown in dashed curve and dotted line. 
We use the same regions where the bulg e and disk clearly dominate the observed profile as 
adopted in the paper of iBorosonI (Il98ll ) to derive the decompositions iteratively: between 
3.5" and 23.0" for the bulge and between 69.0" and 230.2" for the disk. The decomposition 
procedure gives the central surface brightness (/x^) of 19.64 ± 0.19 mag arcsec"^, the effective 
radius (r^) of 0.41± 0.04 kpc (0.16 ± 0.02 arcmin), the Sersic index (n^) of 0.82 ± 0.08 for the 
Sersic law of the bulge, the disc central surface brightness /i^ = 21.19 ± 0.11 mag arcsec"^ 
and scale length h = 3.25 ± 0.20 kpc (1.30 ± 0.08 arcmin) for the expon ential law of th e 
disk. We ob tained a much s n aaller effective radius of the bulge than that of iBorosonI (jl98ll ). 
because the Ide Vaucouleurd (119481 ) law can not fit the bulge component very well when we 
check Figure 6 of their paper . Although larg e discrepancy for the bulge, the scale length of 
the disk approximate those of lBorosonI (jl98ll ) corresponding to th eir adopted distance of 12.2 
Mpc. Actually, m any exceptions when fit ting the bulge with the Ide Vaucouleurd (119481 ) law 
were mentioned by lMcDonald et al.l (120091 ) (and references therein), and it is m ore acceptable 
that the projected three-dimensional bulge profile is fitted by the generalized ISersid (119681 ) 
law. 



Although blue bands are used to derive the structural parameters of different compo- 
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nents traditionally, their light profiles are substantially affected by the recent star formation 
and the extinction of the gas and dust. To the contrary, the NIR luminosity is a good tracer 
of the stellar mass and maps the distribution of the old stellar populations in galaxies. It 
is insensitive to the luminosity of young massive stars and unaffected by the extinction of 
the gas and dust for its ignorable extinction coefficient. Therefore, we also decompose the 
surface brightness profile of the NIR Kg band into two components as shown in the right 
panel of Figure |3l The central surface brightness, the effective radius, and the Sersic index 
of the bulge are 16.87 ± 0.18 mag arcsec'^, 0.52 ± 0.04 kpc (0.21 ± 0.02 arcmin), and 1.31 
± 0.11, respectively. The disc central surface brightness and scale length are 19.35 ± 0.10 
mag arcsec"^ and 2.51 ± 0.12 kpc (1.00 ± 0.05 arcmin), respectively. The effective radius, 
Ser sic index, arid disc scale length are consistent with the structural parani eters derived 
by iGanda et al.l ( 120091 ) for the HST H band photometry ( IGanda et al.ll2009[ ). which give 
re = 0.21 d= 0.002 arcmin, rib = 1-23 ± 0.01, and h = 1.18 ± 0.02 arcmin. From the fitting 
parameters, we find that the disk profile decrease more dramatically in the NIR band than 
in blue optical band, since there are a number of young stellar populations located in the 
disk and they raise relatively more luminosity in the blue bands. The effective radius of the 
Kg band is larger than that of the d band, giving a more reliable span of the stellar mass 
of the bulge. We regard the region of R < 0.21' as the region of the bulge, the region of 
R > 0.5' where the profile starts to deviate from the exponential law as the region of the 
disk, and the region between these two RrGclS clS db transition zone. 



5.2. Age Distribution 



Figure E] displays the age distributions of NGC 628 in 2D (left panel) and histogram 
(right panel). In the left panel of Figure El we find that age is oldest near the galactic 
core and becomes younger with the increasing radial distance (the age ranges from about 10 
Gyr to 2.0 Gyr). Two apparent spiral arms in the disk, where considerable Hll regions are 
concentrated, are much younger (~ 1.3 Gyr) than any other components. They extend from 
the center to the outer region of this galaxy. The inter-arm areas are filled with relatively 
older stellar populations. 

In order to derive the statistical properties of the stellar populations in more physical 
meaning, the near-infrared luminosity is used as the mass weight for each pixel in the galaxy 
image. Since the near-infrared luminosity is insensitive to young massive stars, and reddening 
in this wavelength range can be also ignored, the nea r - infrared bands, espec ially the Kg band, 
are ideal tracers of the stellar mass ( ICole et al.ll200ll : iKochanek et al.ll200ll ). In the following 
of this paper, we will derive various statistical features (e.g., average values, histograms, and 
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radial profiles) of age, metallicity, and reddening weighted by the Ks band fiux. 

The average age of the whole galaxy, bulge and disk is about 4.9 Gyr, 7.5 Gyr and 4.4 
Gyr, respectively. Half of the total stellar mass is younger than about 3.5 Gyr as seen in the 
right panel of Figure |5l In this figure, we can also see that there are quite a few extremely 
young stellar populations with age less than 1 Gyr. A majority of them are located close to 
the double spiral arms and some Hll regions as confirmed by their positions on the age map. 



Sanchez et al.l ( 1201 ll ) presented a wide-field IFS survey on NGC 628 and obtained 2D 
spectra in the FOV of about 6 arcmin in diameter with spatial resolution of about 2.7 arcsec. 
A dithering mode was adopted for the central pointing to increase the spatial resolution. 
They derived the global age of about 8.2 Gyr by both fitting the integrated spectrum of the 
whole observed galaxy to linearly combined SSPs and averaging the radial age distribution. 
The discrepancy of the average age with respect to our estimates mainly originates from the 
different size of the o bserved galactic area, which in our study is more than 2 times that of 



Sanchez et al.l ( 1201 ll ). So considerable rather young stellar populations in the outer region 



of the disk are involved in the calculation of our average age. 



5.3. Two Distinct Disc Components 

The radial age profile in the left panel of Figure |6] shows that there are two distinct 
components in the disk, where we consider the radius range of 0.5 - 1.0' as the inner region 
of the disk and the range from 1.0' to 2.2' as the outer region of the disk (separated at the 



radius of about 1.0'). This radial age profile is very similar to that of ISanchez et al.l (12011! ) 
within the radius they analyzed regardless of the system uncertainties. The inner region of 
the disk is older and its age gradient is much steeper than that the outer region. The mean 
ages of these two components are about 7.2 Gyr and 2.9 Gyr, respectively, and the slopes are 
11.9 ± 1.0 Gyr arcmin"^ (4.8 Gyr kpc~^) for the inner region in the solid line and 1.2 ± 0.1 
Gyr arcmin"^ (0.5 Gyr kpc~^) for the outer region in the dashed line of the left panel of 
Figure El 

Such distinct disc components in the age distribution should appear in the color profile 
which is related to the properties of the underlying stellar population, such as age and 
metallicity. We calculated the radial profile of the BATC d — g color as plotted in the right 
panel of Figure El These two bands are selected according to their effective wavelengths 
approximating the broad bands of B and V . The variation of this color profile is similar to 
that of the age profile: a fiattening or even inverse tendency within 0.5' and two different 
gradients in the region of i? > 0.5'. 
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Although the inner and outer regions are both in the exponential disk, the age of the 
inner region is much older than that of the outer one. The steeper age gradient of the inner 
region of the disk indicates this smaller component should experience a very long evolution 
history in an age span of about 7 Gyr. The outer region of the disk with much shallower 
gradient might be formed 2-3 Gyr ago within a time interval of about 1 Gyr, which implies 
that NGC 628 might endure some encounter or accretion events and large number of gas fell 
in to produce new stars and generate the large-scale outer region of the disk in a very short 
time (see Sectioning])- 



5.4. Young Circumnuclear Ring Structure 



In the left panel of Figure [6l contrary to the disc components, an inverse age gradient 
resides in the bulge and the transition zone (also discovered in the d — g color profile), while 
the oldest age up to 10 Gyr comes up at the galactocentric radius of about 0.5' (1.3 kpc) 
instead of the galactic center. Relatively younger structure near the galactic core indicates 
that some processes like gas inflow directly through the galactic disk to the bulge might occur 
in the history and trigger the star formation. This may be also associated to the possible 
existence of a weak barr ed pot e ntial which induces the circumnuclear star formation as 
presented in the paper of ISeigarl (120021 ). 



Natali et al.l (119921 ) showed the broad band UBVRI surface photometry of NGC 628 
and found an obvious blue drop in their B — V profile for i? < 1.5 kpc, which can be predicted 



by mo dels due to infall of gas in the bulge. Sub-mm CO (l- O'l observation oflWakker fc Adler 
(119951 ) and infrared 2.3 /im CO absorption spectroscopy of I James fc Seigarl (119991) displayed 



the ex istence of a circumnuclear star forming ring in the center of the galaxy. iCornett et al. 
(I1994J ) used the UV surface photometry and discovered that the nuclear region has the 
overall morphological characteristics of spiral arm materials resembling those of M33 and 
might endure a significant star formation in the past few Gyrs predicted by the UV-optical 
colors. 



In the paper of ISanchez et al.l (120111 ). in addition to the decreasing age gradient in the 
inner region and possible shallow gradient in the outer region of the disk, an inverse age 
gradient was also obtained in a circumnuclear ring at about 25", which is very close to the 
position of our result. This kind of ring structure was also detected by iGanda et al.l (120061 ) 
in their H/3 and [O ill] distributions, obtained by the integral field spectroscopy with the 
SAURON Integral Field Unit pointing to the galactic core (about 33"x41"). 
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5.5. Metallicity Distribution 



The metallicity Z, derived by SEP fit, is conv erted to [Fe/H] using the standard chemical 
composition of the Sun fiGrevesse &: Sauvallll998l ). Related histogram and radial distribution 
are also calculated. In the left panel of Figure [71 we can see that the metallicity is lower 
near the core than in the outer regions and the abundance close to spiral arms is much richer 
than the inter-arms due to their poles of the frequently born place of massive stars. 

Average [Fe/H]s of the whole galaxy, bulge, disk, and the inner and outer region of 
the disk are about -0.77, -1.04, - 0.78, -0.93, and -0.71 dex, respectively. The derived global 
metallicity is lower than that of ISanchez et al.l (120111 ). whose spectra mostly lie along the 
spiral arms and in the bulge. Much r i cher a bundances of Hll regions in the spirals make 
the global metallicity of ISanchez et al.l (l201l[ ) larger than ours. In the histogram of [Fe/H] 
(middle panel of Figure [7]), two components dominate this spiral galaxy: the older stellar 
population with poorer abundance and younger population with richer metallicity. The 
poor stellar population mostly lies within the galactic core and inter-arms, while the rich 
one mainly resides close to the spiral arms and the outer region of the disk. The distributions 
of these two populations take on a bimodal profile. We fit the overall distribution by a double 
Gaussian function which has two single Gaussian functions as its items. Parameters of these 
two Gaussian components are summarized in Table IH The two components and the bimodal 
curve are respectively plotted in solid and dashed lines in the middle panel of Figure [71 



As presented by Hll observations of iTalentI ( ll983l ) and lBelley fc RoyI (Il992[ ). an obvious 
radial 0/H abundance gradient exists in the disk of NGC 628, which is 0.21 dex arcmin"^ and 
0.17 ± 0.004 dex arcmin"^ , respectively corresponding to their adopted distances. Equally, 
a weak decreasing [Fe/H] gradient located in the outer region of the disk is shown in the 
right panel of Figure [71 [Sanchez et al.l (120111 ) also obtained the radial stellar metallicity 
distribution, the profile of which is very similar to that of ours. The slope of this weak 
gradient is about 0.11±0.03 dex arcmin"^ (0.05 dex kpc~^). Comparing those gradients with 
our result of 0.11 dex arcmin"^, we discover that the stellar metallicity gradient is somewhat 
lower than that of the gas abundance. It is notable that there is a reverse gradient of the 
inner region of the disk with a slope of 0.48 ± 0.08 dex arcmin"^ (0.19 dex kpc~^). Such 
kind of abundance distribution in the disk provides fundamental constraints on the chemical 
evolution of NGC 628. 



-17- 



5.6. Intrinsic Reddening Distribution 



Both dust and gas absorb and scatter star light and make the galaxy redder. Extinctions 
are much larger near the UV bands than in the IR bands. Accurate photometry from the UV 
to infrared can derive reliable intrinsic reddening {E{B — V)) of NGC 628 as shown in the left 
panel of Figure [SI In this map, regions close to the spiral arms and the bulge hold relatively 
larger reddenings. The mean reddening of those regions close to spiral arms is about 0.38 
and show no r adial gradient, which is consistent with the observational measurements in 
Belley fc RoyI (119921 ). Some inter-arm regions and parts of the outmost disk appears to be 
less reddened by dust and gas. 

In order to confirm the reliability of the reddening value, dust emission of the Spitzer 
8.0 fim is compared with the reddening distribution. We construct the intensity map from 
the dust in Spi tzer S.Oum ( r ight p anel of Figure [8]) dominated by the PAH emissions using 
the method of iHelou et al.l (120041 ). who subtracted the background stellar fluxes from the 
obse rved fluxes. I n the determination of the stellar fluxes in longer wavelengths than 3.6 
(120041) extrapolated the 3.6 /xm flux using the stellar population model 



Helou et al. 



of StarBurst99 (ILeitherer et al.l Il999l ) to obtain the scale factor of 0.232 in 8.0 /im (i.e., 
Dg.o = h.o — 0.232/3.6, where Dg.o is the intensity of the dust emission and /s.e and Is.o are 
the observed intensities in the Spitzer 3.6 and 8.0 /xm, respectively). Superposed contours 
in the right map of Figure [H] show the distribution of 2.6 mm CO (1-0) emission intensity 
from the Berkeley- IUi nois-Maryland Asso ciation Survey of Nearby Galaxies (BIMA SONG) 
(IRegan et al.ll200ll ). iRegan et al.l (120061 ) concluded that the CO and 8 fim emission mor- 
phologies are similar as also presented in this map and 8 fim PAH surface brightness can be 
used as a possible tracer of the interstellar medium. Clearly strong PAH emissions prefer- 
entially distribute along the spiral arms where a number of Hll regions radiate UV photons 
to photodissociate the PAHs around the molecule clouds and generate IR emissions. Com- 
pared with the reddening map, we discover that many of the regions with large reddening 
values in the galaxy show luminance peaks in the 8 /im PAH image, implying the existence 
of abundant gas and dust. Some other regions with large reddenings might be caused by the 
absorption of dust lanes as checked in the optical observations with respect to the positions 



of pri mary dust lanes and reddening feathers, which were investigated by iLa Vigne et al. 
J2006h . 



Histogram in the left panel of Figure M shows that E{B — V) varies from 0.1 to 0.6 
and the mean reddening is about 0.31. Average reddenings of the bulge, the disk, and the 
inner and outer regions of the disk are about 0.41, 0.31, 0.32, and 0.31 mag, respectively. 
The radial distribution in the right panel of Figure [9] presents a very weak gradient of the 
intrinsic reddening with a slope of 0.015 ± 0.002 mag arcmin"^. Several bumps in this plot 
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are produced by the large reddening values near the spiral arms. We can see in the radial 
profile that the bulge as well as the transition zone has larger reddening values than the 
outside area, supporting the young circumnuclear ring structure holding abundant gas and 
dust. That is, a considerable volume of dust and gas should be depos ited in the c ore o f 
NGC 628, giving a rise of star formation in its recent history. In contrast, iKong et al.l (120001 ) 
studied another nearby galaxy (M81) of Sab type using a similar method. Their results show 
that M81 has a much lower reddening range of 0.08 - 0.15 mag in the bulge than ours and 
the disk reddening (~ 0.2 mag) is also somewhat smaller, implying that NGC 628 might be 
a gas-rich galaxy. 



6. Discussions 



6.1. Pseudobulge and Secular Evolution 

The Sersic indices in both optical and NIR bands as we have derived from the bright- 
ness profile fitting are close to 1, which represents the exponential law. It has been in- 
creasingly discovered in recent years that bulges with the surface brightness profiles close to 
the exponential law are sig nificantly different from those with larger Sersic indices resem- 



bling the elliptical galax i es (Kormenc 



2008; Ganda et al. 2009; Fisher et al. 



yl 1993 : Kormendy &: Kennicutt 2004; Fisher fc Drory 



20091 ). Now, two types of bulges are typically called: 



classical bulges that are dynamically hot and featureless and pseudobulges that retain the 
memory of the disk origin. The pseudobulge s are v ery common in late-type galaxies (75% 
of 77 Sd-Sm galaxies as given in iBoker et al.l (|2002[ )). They have flatter shapes, kinematics 



dominated by rotation and corresponding smaller velocity dispersion, active star formation, 
nuclear bars, nuclear rings, and/or nuclear spirals and nearly exponential brightness profiles 
[nb ~ 1), while the classical bulges are dominated by r andom motions, contain old stellar 



popu l ations, and are mo re similar to the E-type galaxies ( Fisher fc Drorv 2008[ 



Fisher et al 



Ganda et al. 



J2009h i dentified 



20091 : iFisher et al.ll2009l . and references therein). Specifically, 
pseudobulges as those bulges containing nuclear bars, nuclear spirals, and/or nuclear rings 
and having the Sersic index less than 2. 

The Sersic indices of NGC 628 are 0.83 for the blue band and 1.31 for the NIR band, 
both of which approxima te the exponential l aw. Nuclear spiral arms are clearly seen in the 
UV band of this galaxy (ICornett et al.lll994l ). Moreover, the young ring structure detected 
in our age map and cental concentration of rich gas shown in the reddening map declare 
recent young star formation in the galactic center. Therefore, the bulge of NGC 628 belongs 
to that so-called disk-like pseudobulge. 
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Unlike the classical bulges, which are typically merger-built, ps eudobulges have the op- 
portunity to be grown v ia the internal secular evolution of the disk ( iKormendy &: Kennicutt 
2004 : iFisher et al.ll2009[ ). Disk galaxies can evolve by the rearrangement of mass and angular 
momentum driven by the non-axisymmetries of bars, ovals and/or spiral structures, which 
actuate the gas infall towards the galactic center, t rigger the star formation , and build up the 
central mass concentrations of the pseudobulges (IKormendy &: KennicuttI l2004j ) . In barred 
or oval disk galaxies, gravitational torques and/or shocks near their potential minimum cause 
the gas inflow and outflow, which can bring together the disk gas, hence trigger the star for- 
mation and make outer rings, inner rings , and central mass concentrations (pseudobulges) 
( iRegan &: Teubenll2003l : iFisher et al.ll2009l ). In normal spiral galaxies without bars and ovals, 
spiral structures are maintained by the density wave that propagate through the disk. Shocks 
are formed as gas approaches and leaves the arms, making the gas lose the energy and sink 
into the center, and form the disk-like bulge (IKormendy &: KennicuttI 12004 ) . NGC 628 is a 
late-type spiral galaxy without evident bars and ovals, so it is possible that its pseudobulge 
was formed by the secul ar evolution of the disk driven by the non-axisymmetric potential 
of spiral arms. Actually, IFisher et al.l (120091 ) concluded that the median pseudobulges could 
have grown the current stellar mass at their present-day SFRs in 8 Gyr and their results are 
consistent with a scenario in which bulge growth via internal secular star formation. The 
age of the ps eudobulge in NGC 628 as we obtained previously is about 7.5 Gyr, similar to 
the results of IFisher et al.l ( 120091 ) and hence supporting the secular evolution. 



Besides secular evolution, possible effects that build the pseudobulges include extremely 
gas-rich accretion events, distant gravitational encounters, an d gravitational interactions in a 
cluster ( IKormendy fc KennicuttI 12004 : IFisher fc Droryll2008l ). NGC 628 is an isolated galax- 
ies and no la rge galaxies are found with in close distance which could generate possible tidal 
interactions (IKamphuis fc Briggsl Il992l ). Thus, gravitational encounters and interactions 
sh ould be unimportant f or the formation of the pseudobulge in this galaxy. Hi observations 
by iKamphuis fc Briggsl (119921 ) show two symmetrical high velocity complexes (HVCs) and 
an extended tail to the southwest. Gas accretion events from an external Hl-rich object may 
form the oute r tail, and accretion c an als o create extreme velocity deviations making those 
HVCs visible. IKamphuis fc Briggsl (Il992[ ) argued that we might be witnessing the flattening 
of the disk after accreting nearby companions long time ago according to the orbit of the tail, 
disc warp sustained by a steady gas infall, and complicated velocity distortion in the outer 
disk. However, it is extraordinarily uncertain whether this accretion-induced formation of 
pseudobulges can retain the disk-like properties. They also did not exclude the primordial 
origin that regards the observed gas distribution as the left-over of the initial formation of 
the galaxy. Nevertheless, the primordial origin and external accretion of gas disk in NGC 
628 do not prevent this galaxy from forming the pseudobulge by the processes in the secular 
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evolution and the secular evolution can still have an effect. 



6.2. Nucleus 



Both the age map (Figure [5]) and the reddening map (Figure [HI) show the nucleus with 
young age (about 3.5 Gyr) and rich gas (about 0.5 mag in E{B — V)). Decomposition of 
the NIR H band in extremely high resolution of the HST image also presents that the radial 
brightness profile e xceeds the Sersic law of the pseudobulge in the innermost of the galaxy 
( IGanda et al.ll2009[ ). This kind of compact nuclei h osts nuclear star clusters and they are 
very common in late-type galaxies ( iBoker et al.ll2002[ ). The nuclei are obviously distinguished 
from the ambient (pseudo)bulges and disk in the sense of a much smaller effective radius 
and higher effective surface brightness. Stellar populations of nuclei tend to have blue colors 
i mplying young age for the stars which contribute most of the light (typical example of M33 



m 



Long et al.l ( 12002| )). The nucleus and the pseudobulge can both exist in the same galaxy 



but they show a great number of differences, making the nucleus harder to be understood. 
Maybe these tiny and dense nuclei we re formed via dynamical frictio n driving the clusters 
sinking to the center as mentioned in iKormendy &: KennicuttI (120041 ). However, the origin 
and evolution of nuclei are quite complicated and its discussion is far beyond the ability with 
current results of our study. 



6.3. Scenarios of Disk Evolution 



In our results, we discover that the disk of NGC 628 can be divided into two parts: 
the inner region and the outer region. These two parts have different properties of stel- 
lar populations as far as age and chemical abundance be concerned. Disk dynamics points 
that the inner part of the disk shrinks and the outer part expands due to angular momen- 
tum transport cause d by differential rotation and non-axisymmetric self-gr avitating mode, 
such as spiral arms ( iLynden-Belll Il979l : iLin fc Pringld 119871 : iTremaind Il989l ) . This redistri- 
bution of mass and angular momentum results from the minimization of the total energy 
with the total angular momentum conservation. The gas contraction in the inner part of 
the disk could more easily trigger the star formation than in the outer part, and hence 
evolve more quickly. Merely by this kind of internal secular evolution, the disk can con- 
struct two different regions with a large age discrepancy. Certainly, one should confirm 
whether this process is fast enough to form the structures within the Hubble time. Non- 
axisymmetry of the pote ntial caused by the spiral arni s can provide the engine for rapid 
evolution as discussed in iKormendy fc KennicuttI ( 120041 ). On the other hand. Hi observa- 
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tions of IShostak &: van der KruitI (Il984l ) and iKamphuis fc Brigga (Il992l ) display that most 
of the gas concentrates on the disk and an extensive tail lies to the southwest. Gas accretion 
events from an Hl-rich object with Mhi ~ 9 x IO^Mq m ight occur long time ago and they 
affect only the outer disk loosely bounded by gravitation (IKamphuis fc Briggslll992l ). There- 
fore, the accreted gas could rule the following star formation activities in the outer region, 
forming a low stellar-density and young outer region of the disk. 



7. CONCLUSION 

Recently, evolutionary stellar population synthesis models become more and more pop- 
ular in studying the properties of the stellar populations in order to understand the star 
formation and evolution histories of galaxies. In this paper, we adopt one of the EPS mod- 
els, PEGASE, to analyze the age, metallicity and reddening distributions of a nearby face-on 
spiral galaxy (NGC 628) based on multi-band photometric data ranging from UV to IR. On- 
line archival data from GALEX, XMM OM, 2MASS, Bok telescope and Spitzer provide the 
ultraviolet and near-infrared photometry of this galaxy. The optical observations are ob- 
tained in the BATC survey with 15 intermediate bands ranging from 3000 A to about 1 /im. 
By comparing the photometric SED (totally 23 bands) with those of model ones calculated 
from the PEGASE synthetical spectra, we derive the properties of stellar population for each 
part of the galaxy (area of 1.7" in the BATC pixel size) and hence the spatially resolved age, 
metallicity and reddening maps. 

Structure parameters are calculate d in optical and NIR bands by fitting the surface 



brightness profiles with the ISersid (119681 ) law for the bulge and exponential law for the disk, 
in order to discuss the properties of different components in NGC 628. The bulge size is 0.52 
kpc and disc scale length is 2.51 kpc. The mean stellar age of the galaxy is about 4.9 Gyr, 
so young that a great number of stars might form in the recent history of 2 - 3 Gyr. The 
age of the bulge is about 7.5 Gyr, presenting the oldest stellar population. The radial profile 
shows the age becomes younger and younger from the center to the outer disk. We confirm 
the presence of a young circumnuclear ring structure in the bulge and the transition zone 
within the galactocentric distance of about 0.5'. This young structure behaves to have flat 
and even inverse age gradient, showing active star formation. The disk can be divided into 
two parts as presented in the radial age profile: the old inner region and young outer region, 
which have distinct age gradients of 11.9 Gyr arcmin"^ and 1.2 Gyr arcmin"^, respectively. 
The mean ages of these two regions are about 7.2 Gyr and 2.9 Gyr. Both the circumnuclear 
ring structure and two different regions of the disk can be clearly identified in the BATC 
d — g color profile. 
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Abundance map of NGC 628 in [Fe/H] reveals that the bulge is consisted of old stellar 
populations with low metallicity ([Fe/H] ~ -1.04 dex) and the disk is relatively richer ([Fc/H] 
~ -0.78 dex). A very weak abundance gradient with slope of 0.11 dex arcmin"^ (0.05 dex 
kpc~^) for the outer region of the disk is obtained in the radial profile, which is smaller than 
the gas abundance gradient obtained by the measurements of Hll regions. Much larger [Fe/H] 
values in the metallicity map are discovered neighboring to Hll regions, where numerous 
young massive stars have been being formed. Reddening map shows obvious spiral-like lanes 
extending from the galactic core to the outer disk. Comparing with the IRAC 8 yum dust 
emission image, we find that many parts of the reddening distribution correlate the PAH 
emissions. The global average reddening in E{B — V) is about 0.31. Reddening values of the 
bulge and disk are about 0.37 and 0.31, respectively. A tiny gradient of the radial reddening 
distribution is found, whose slope is 0.015 mag arcmin"^. 

At the end of this paper, the formation and evolution of different components of NGC 
628 are discussed based on the results. As very common in late-type galaxies, the bulge of 
this galaxy is a pseudobulge, which has the Sersic index close to the exponential law, young 
star formations and central spiral structures as seen in the UV band. The spiral structures 
generate the non-axisymmetric potential, make the gas fall into the center and form the 
central gas concentration which triggers the star formation and grows the pseudobulge. 
All these processes are secular and our results about the pseudobulge support its growth 
by the secular evolution. However, Hi gas distribution in the disk does not rule out the 
primordial origin and gas accretions affecting the formation of the pseudobulge. We also 
discover that a nucleus hosting star clusters can be obviously noticed, which has very young 
stellar populations and large reddening. Since NGC 628 is an isolated galaxy, gravitational 
interactions with other external galaxies are excluded in both the bulge and disk evolution. 
The secular evolution of the disk can form two distinct regions with a great age discrepancy, 
although gas accretion if existing, could also affect the outer region. 
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Table 1: Filter information in the BATC photometric system and statistics of observations 



No. 


Filter 


Acff" (A) 


Bandwidth(A) 


Exp^ (s) 


FWHM'^ (") 




rills'^ 


01 


a 


3360 


222 


13200 


4.61 


2 


0.041 


02 


b 


3890 


187 


8400 


4.45 


3 


0.016 


03 


c 


4210 


185 


7200 


4.01 


6 


0.008 


04 


d 


4550 


222 


11400 


3.43 


5 


0.015 


05 


e 


4920 


225 


13800 


3.73 


3 


0.009 


06 


f 


5270 


211 


11100 


4.04 


3 


0.006 


07 


9 


5795 


176 


6000 


4.44 


2 


0.002 


08 


h 


6075 


190 


6000 


3.48 


4 


0.004 


09 


i 


6660 


312 


3900 


4.28 


5 


0.005 


10 


j 


7050 


121 


6600 


4.24 


7 


0.005 


11 


k 


7490 


125 


7500 


5.07 


5 


0.009 


12 


m 


8020 


179 


7800 


5.25 


1 


0.021 


13 


n 


8480 


152 


10800 


4.44 


5 


0.005 


14 


o 


9190 


194 


15300 


4.57 


8 


0.009 


15 


P 


9745 


188 


19900 


4.89 


1 


0.020 



''Effective wavelength, of eacli filter 

''Total exposure time of the combined image 

^Full width at half maximum (FWHM) of the combined image 

''image number for calibration 

"^Calibration error in magnitude as explained in Section l3.ll 
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Fig. 1. — Images of NGC 628 in three BATC bands centered at the wavelengths of 3890 
(left), 6660 (middle), and 9190 A (right). The image scale is 361x361 (1.71'7pixel). 



- 30 - 



Table 2: Multi-band observations from different telescopes for NGC 628 



Name'' 


Filter 


(A) 


Bandwidth 

(A) 


FOV* 


scale° 
(") 


Exp* 


FWHMS 
(") 


Reference 


GALEX 


FUV 


1516 


268 


1.28° 


1.5 


1636.05 


4.3 


1 


NUV 


2267 


732 


1.24° 


1.5 


1636.05 


5.3 




UVWl 


2905 


620 


17' 


0.95 


~1000.0 


2.0 




XMM-OM 


UVM2 


2298 


439 


17' 


0.95 


~1000.0 


1.8 


2 




UVW2 


2070 


500 


17' 


0.95 


-1000.0 


2.0 




BATC 


a - p 


3000 - 9900 


120 - 310 


58' 


1.7 


1.1 - 5h 


~4.3 


3 




J 


1.235 fim 


0.162 fim 


0.39° 


1 


7.8 


~2.5 




2MASS 


H 


1.662 fim 


0.251 Atm 


0.39° 


1 


7.8 


~2.5 


4 




Ks 


2.159 Atm 


0.262 Atm 


0.39° 


1 


7.8 


~2.5 




Bok 


Ks 


2.16 /im 


0.262 Aim 


8.5' 


0.24 


~1500 


1.5 


5 




IRACl 


3.550 fim 


0.75 Aim 


5.2' 


0.75 


1500.8 


1.66 




Spitzer 


IRAC2 


4.493 Aim 


1.01 Aim 


5.2' 


0.75 


1500.8 


1.72 




IRAC3 


5.731 fim 


1.42 Aim 


5.2' 


0.75 


1500.8 


1.88 


6 




IRAC4 


7.872 /im 


2.93 Aim 


5.2' 


0.75 


1500.8 


1.98 





"Name of the telescope or survey 
^Name of the filter 

'^Effective wavelength in angstrom if not specified 
"^Field of view of the CCD camera 
•^Mosaic pixel scale 

•'^ Total exposure time in unit of seconds if not specified 
^Full width at half maximum of the mosaic 



Note. — Referenc es. ( 1) iMorrissev et al. 
Skrutskie et"all |2006l ): (5) iKnapen et al.l (|20q4 )" 



; (6) iFazio et all (|2004D 



( 200/1): (2) iKuntz et al l (120081) : (3) IZhou et al.l (|200lh : (4) 
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R.A.(deg) 

24.30 24.25 24.20 24.15 24.10 24.05 




* * 



1 5°42 




oi'jy'^o" 



36^50= 



'lO" 

R.A.(J2000) 



Fig. 2. — True color map of NGC 628 synthesized with the 8 /im band (red), Ha narrow 
band (6563 A, green) and NUV band (2271 A, blue) images. These bands are chosen to 
present the young stellar population and dust emission. Red and green solid circles are stars 
and Hll regi ons to be re r aoved . Only those regions with radius larger than 5 arcsec in the 
catalogue of iFathi et al.l ( 120071 ) are displayed in this figure. Black squares are the selected 
pixels for illustrating the model fitting goodness in section 14.31 The solid circle around the 
center, the dashed circle, and the outmost dash-dotted circle represent the bulge together 
with the transition zone within R < 0.5', the outer boundary of the inner region of the disk 
at 1 arcmin, and the outer boundary of the outer disk at about 2.2' where the valid SEDs 
are still available (see Section [573]) . respectively. 
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Bulge effective radius; 0.41 kpc 
Disc scale length: 3.27 kpc 
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Fig. 3. — Left: radial surface brightness profile of the BATC d band as plotted in filled 
squares with error bars. Right: radial surface brightness profile of the BATC Kg band as 
plotted in fill ed squ a res w ith error bars. The dashed curve and dotted line in both panels 
are the fitted ISersid ( 1l968l ) law for the bulge and exponential law for the disk, respectively. 
The solid curve is the sum of these two fitted components. 



Table 3: Observed SEDs of several sampled pixels and corresponding fitted parameters 



No 



Y ci(J2000) i5(J2000) NUV UVM2 UVWl a 



f 



rn n 



H Ks IRAQI IRAC2 Age E(B-V) Z [Fc/H] 



1 182 189 24.1739 15 

2 188 205 24.1766 15 

3 244 218 24.2041 15 

4 196 155 24.1813 15 

5 212 128 24.1895 15 

6 134 199 24.1501 15 



7868 22.76 

0.18 
7945 23.48 

0.27 
8014 24.84 

0.44 
7708 23.52 

0.28 
7582 23.82 

0.29 
7909 22.37 

0.16 



22.89 

0.18 
23.90 

0.37 
99.99 

9.99 
23.86 

0.28 
23.62 

0.31 
22.60 

0.18 



21.67 

0.12 
22.65 

0.19 
23.65 

0.30 
22.65 

0.18 
23.07 

0.22 
21.81 

0.12 



20.92 

0.09 
21.91 

0.15 
23.06 

0.33 
21.98 

0.16 
22.56 

0.23 
21.51 

0.11 



20.19 19.78 
0.06 0.05 
21.32 20.98 
0.11 0.09 
22.41 22.28 
0.20 0.21 
21.32 21.04 
0.11 0.09 
21.87 21.60 
0.15 0.13 
21.07 20.98 
0.09 0.09 



19.46 19.15 

0.05 0.04 
20.66 20.38 

0.08 0.07 
21.96 21.79 

0.17 0.14 
20.74 20.49 

0.08 0.07 
21.45 21.18 

0.12 0.10 
20.80 20.64 

0.08 0.08 



18.91 
0.03 
20.13 
0.06 
21.52 
0.12 
20.27 
0.07 
20.97 
0.09 
20.49 
0.07 



18.65 
0.03 
19.94 
0.06 
21.34 
0.11 
20.08 
0.06 
20.86 
0.09 
20.39 
0.07 



18.61 18.42 

0.03 0.03 
19.90 19.72 

0.06 0.05 
21.37 21.25 

0.12 0.12 
20.05 19.86 

0.06 0.06 
20.82 20.67 

0.09 0.08 
20.39 20.20 

0.07 0.07 



18.30 18.12 
0.03 0.02 
19.62 19.48 
0.05 0.05 
21.16 21.04 
0.12 0.12 
19.77 19.64 
0.05 0.05 
20.60 20.49 
0.09 0.09 
20.19 20.07 
0.07 0.07 



18.05 

0.03 
19.40 

0.07 
20.89 

0.28 
19.57 

0.08 
20.43 

0.15 
20.05 

0.11 



17.98 17.82 

0.03 0.03 
19.36 19.20 

0.06 0.07 
20.73 20.61 

0.16 0.18 
19.49 19.35 

0.07 0.07 
20.36 20.30 

0.13 0.16 
19.91 19.77 

0.09 0.10 



17.73 17.54 

0.03 0.02 
19.14 19.09 

0.07 0.06 
20.70 20.28 

0.29 0.17 
19.27 19.28 

0.08 0.08 
20.22 19.98 

0.17 0.12 
19.81 19.43 

0.13 0.08 



17.33 
0.02 
18.68 
0.06 
20.31 
0.26 
18.92 
0.07 
20.01 
0.14 
19.51 
0.10 



17.52 
0.03 
18.95 
0.08 
99.99 
9.99 
19.07 
0.09 
20.08 
0.19 
19.58 
0.14 



18 .32 
0.03 

19 .73 
0.05 

21 .27 
0.10 

19 .83 
0.05 

20 .76 
0.08 

20 .16 
0.06 



18.83 
0.03 
20.24 
0.06 
21.66 
0.12 
20.29 
0.06 
21.24 
0.10 
20.59 
0.08 



13.2 
0.6 
8.5 
0.9 
1.1 
0.3 
3.2 
0.4 
2.6 
0.5 
1.3 
0.2 



0.34 
0.01 
0.30 
0.01 
0.43 
0.03 
0.37 
0.01 
0.27 
0.02 
0.18 
0.02 



0.0012 
0.0001 
0.0016 
0.0003 
0.0010 
0.0005 
0.0016 
0.0003 
0.0032 
0.0010 
0.0302 
0.0040 



-1.18 
0.04 
-1.09 
0.08 
-1.23 
0.20 
-1.04 
0.09 
-0.74 
0.19 
0.30 
0.06 



X and Y arc the pixel position in tlie image with the size of 361 X 361. The positions of these pixels are marked in Figure^ Both right ascension a and declination S arc in degrees. The last 
four columns arc the fitted parameters. There are two rows for each pixel: one row contains the observed SEDs and fitted results and the other includes their errors. Note that magnitudes 
of 99.99 and errors of 9.99 arc set for the bands whose SNRs arc three times lower than the sky background standard deviation. 
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Fig. 4. — Combined plots of observed SEDs, matched model SEDs and corresponding spectra 
for six randomly sampled pixels as described in Table [31 The filled symbols are the observed 
AB magnitudes of specified bands and open circles are the matched model magnitudes. 
Spectra are drawn in grey. Vertical bars are the measured errors of the observed magnitudes 
as defined in Equation [2J 
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Fig. 5. — Left: two-dimensional age distribution of NGC 628. Blank regions are masked 
stars and Hll regions. Right: histogram distribution of mass weighted ages. The cumu- 
lative flux in the vertical ordinate is the total apparent luminosity of the Ks bands in 
10~^°ergs s~^cm~^Hz~^arcsec~^ which represents the stellar mass. 
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Fig. 6. — Left: radial age profile. Right: radial profile of the BATC d — g color, 
and dashed lines present two fitted gradients in the inner region of the disk (0.5 - 
the outer region of the disk (1.0 - 2.2 '), respectively. 
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Fig. 7. — Left: two-dimensional [Fe/H] distribution of NGC 628. Blank regions are masked 
stars and Hll regions. Middle: histogram distribution of mass weighted metallicity. The 
cumulative flux in the vertical ordinate is the total apparent luminosity of the Ks bands 
in 10~^°ergs s~^cm~^Hz~^arcsec~^. Two components of poor and rich metallicity are fitted 
by a double Gaussian function (bimodal), which is plotted in the dashed curve. The solid 
curves are two single Gaussian items of the bimodal. Right: radial [Fe/H] distribution. The 
solid line and dashed lines show the [Fe/H] gradients within the range of 0.5 - 1.0' and the 
range of 1.0 - 2.5', respectively. 
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Table 4: Bimodal parameters of the [Fe/H] distribution fitted by a double Gaussian function 



Abundance 




Error^^ 




Error^- 


Ad 


Error^ 


poor 


-1.15 


0.03 


0.23 


0.03 


1.66x10^ 


0.20x10^ 


rich 


-0.40 


0.08 


0.39 


0.07 


1.13x10^ 


0.09x10^ 



"Mathematical expectation of Gaussian function 
''Error (for all parameters) with 95% confidence 
'^Standard deviation of Gaussian function 
"^Constant coefficient of Gaussian function 
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Fig. 8. — Left: two-dimensional intrinsic reddening distribution of NGC 628. Blank regions 
are masked stars and Hll regions. Right: IRAC 8.0 fim stellar flux-subtracted image in the 
same area as the reddening map. Overlapped contours show CO intensities (contour levels: 
1.5, 3, and 5 Jybeam~^kms~^). 
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Fig. 9. — Left: histogram distribution of the overall intrinsic reddening in E{B — V). The 

cumulative flux in the vertical ordinate is the total apparent luminosity of the Kg bands in 
]^Q-30gj.gg s^^cm^^Hz~^arcsec~^. Right: radial distribution of the intrinsic reddening. The 
solid line presents the fitted radial gradient. 



